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Abstract 

In order to model pions of two-flavor QCD we consider a lattice field theory involving two flavors 
of staggered quarks interacting strongly with U(l) gauge fields. For massless quarks, this theory 
has an SUl{2) x SUr(2) x Ua(X) symmetry. By adding a four-fermion term we can break the 
Ua{X) symmetry and thus incorporate the physics of the QCD anomaly. We can also tune the pion 
decay constant F, to be small compared to the lattice cutoff by starting with an extra fictitious 
dimension, thus allowing us to model low energy pion physics in a setting similar to lattice QCD 
from first principles. However, unlike lattice QCD, a major advantage of our model is that we can 
easily design efficient algorithms to compute a variety of quantities in the chiral limit. Here we 
show that the model reproduces the predictions of chiral perturbation theory in the e-regime. 



1 



I. INTRODUCTION 



One of the outstanding problems in lattice QCD is to compute low energy hadronic 
observables, which are dominated by the physics of light quarks, with controlled errors. 
Unfortunately, with current algorithms it is difficult to perform calculations at realistic quark 
masses. In particular, one would like to have light pions of mass 140 Mev, which means 
m 1 ~ Is' w here mi refers to the average up and down quark masses and m s refers to the 
strange quark mass. In practice the lightest pions with current algorithms can reach about 
250-300 MeV with Wilson type fermions [l| and slightly lighter with staggered fermions. 
For example the MILC group has attained — ~ ~ with staggered fermions while — > ~ 
. Given that calculations are usually performed at unphysically large quark 
masses, chiral extrapolations, based on chiral perturbation theory, are used to obtain results 
at realistic quark masses. However, for such an approach to be reliable, we must know the 
range of quark masses over which the chiral expansion used in the extrapolation is valid. 
Unfortunately, this is not yet well understood. 

One way to check if chiral extrapolations are reliable is to verify them in qualitatively 
distinct regimes. Fortunately, while there are many such regimes depending on the values 
of the pion masses and the physical box sizes, the most popular examples are the p-regime 
and the e-regime. Interestingly, in all the regimes the extrapolation formulas depend on the 
same low energy constants that describe the chiral Lagrangian. Hence if the data from a 
lattice calculation can be fit in both the regimes with the same low energy constants one 
would gain more confidence in the extrapolations. As far as we know such comparisons 
have not yet been made in the context of lattice QCD. In fact, although chiral perturbation 
theory is a widely understood and accepted tool, much less effort has gone into describing 
the low energy physics of a fundamental lattice field theory using the chiral Lagrangian. 
Apart from lattice QCD, as far as we know, the only known example where lattice results 
were understood with chiral perturbation theory was in the context of the non-linear sigma 
models and quantum spin systems |6|. However, even in such studies comparison 

between the p-regime and the e-regimes were never made. Further, these studies are some 
what outdated given the advances in computing resources. More such studies and with 
simpler models that resemble QCD closely may teach us more about chiral extrapolations. 
Motivated by this, we introduce and study a lattice field theory model of pions in two-flavor 
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QCD. 

Our model is nothing but strongly coupled lattice QED with two flavors of staggered 
fermions. This model exhibits the same symmetries as two flavor QCD and was recently used 



to study the chira 
in the eighties 



rase transition J7J. The st rong coupling limit was studied considerably 
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171 ] because some of the qualitative physics 



remains despite this limit suffering from the worst lattice artifacts. In fact, even a (7(1) 
gauge theory exhibits confinement and chiral symmetry breaking in the strong coupling 
limit. Notice further that the taste symmetry of continuum staggered fermions is irrelevant 
since it is maximally broken at strong couplings. The relevant symmetry is due to the two 
flavors and resembles QCD closely. Thus, our model allows us to model pion physics from 
a fundamental lattice field theory, very similar to QCD. 

Since the pion is the lightest QCD bound state and plays a crucial role in chiral symmetry 
breaking, even a qualitative description of it from first principles may prove useful. However, 
unless we find a way to fine tune our model, the pion decay constant F (defined here in the 
chiral limit), is naturally close to the cutoff. Thus, the low energy pion physics in our model 
may be described by a chiral Lagrangian that contains information about lattice artifacts. 
In order to circumvent this problem, we define our model in d + 1 dimensions where d = 4 
is the space time dimensions. The extra dimension plays the role of a fictitious temperature 
which allows us to tune F to values much smaller than the cut-off. Thus, we can still explore 
"continuum- like" physics even in the strong coupling limit. 

In addition to its similarity with QCD, another important motivation for studying the 
above model is that the gauge dependent degrees of freedom can be integrated over, which 
significantly simplifies the theory. Further, a new class of algorithms, called the Directed 
Path Algorithm, allows us to study the chiral limit very efficiently Here we extend the 
algorithm to our model and show that the results in the chiral limit agree with predictions 
in the e-regime. A preliminary version of this work can be found in 19]. 
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II. MODEL AND SYMMETRIES 



The Euclidean space action of the model we consider is given by (note that the usual 
factor of | in the fermion kinetic term has been absorbed into the fields): 

d+l 



X [1=1 

where x denotes a lattice site on a d + 1 dimensional hyper-cubic lattice L t x L d . Here 
L t represents a fictitious time direction and will be used to tune F, the non-perturbative 
physical mass scale, to be small compared the lattice cutoff. On the other hand L d will be 
the usual Euclidean space-time box. if> x and ip x are two component Grassmann fields that 
represent the two quark (u, d) flavors of mass m, and 4>n,x is the compact U(l) gauge field 
through which the quarks interact. Here // = 1, 2, d, d + 1 runs over the d + l directions. 
The fj, = 1 direction will denote the fictitious temperature direction, while the remaining 
directions represent Euclidean space-time. The usual staggered fermion phase factors rj^x 
obey the relations: r)f x = T and r)f x = 1 for i = 2,3, d+l. The parameter T controls 
the fictitious temperature and will be used to tune to the continuum limit (F <C 1). The 
coupling c will set the strength of the anomaly. 

When c, m = 0, the action exhibits a global SUl(2) x SUr(2) x Ua(1) symmetry like 
Nf = 2 QCD. To see this it is first useful to note that every lattice site can be classified 
as an even or odd site. Then it is easy to see that the action is invariant under the follow- 
ing transformations: ip — > if} Vlexp(i<f)),ifj — ► exp(i<j))Vnif} , ip e — > ip e V^exp(—i(j)),il) e — ► 
exp(— i(fi)VLip e - Here Vl and Vr are SU{2) matrices and can be parametrized by: exp(i9 ■ a) 
where a\ is a Pauli matrix that acts on the flavor space. At c 7^ 0, Ua{1) is explicitly broken 
and the action is invariant under SUl(2) x SUr(2) x Z2. Thus, the coupling c induces the 
effects of the anomaly. Further at m 7^ 0, it is necessary to set Vl = Vr for the action to 
remain invariant. Thus, with a mass term the chiral symmetry SUl{2) x SUr(2) is explicitly 
broken down to SUy(2). In order to mimic the real world with u,d quarks we need to set 
c 7^ and m/0. Thus, our model has the same chiral symmetry as Nf = 2 QCD. Further, 
based on previous mean field strong coupling calculations 10] , we expect that the symmetry 
breaking pattern is also similar to full QCD. 
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III. MONOMER-DIMER-PION LOOP-INSTANTON REPRESENTATION 



The partition function of our model is equivalent to that of a classical statistical mechanics 
model involving configurations made up of gauge invariant objects such as monomers, dimers, 
pion loops and instantons 20|. We call these MDPI configurations and denote the set of all 
such configurations by /C. Each MDPI configuration is characterized by three site variables 
I (x) , n u (x) , n d (x) and three bond variables 7r"(x), irf(x), ttUx). Here n u (x) is the number 
of uu monomers, rid(x) the number of dd monomers and I(x) the number of instantons 
(or uudd double-monomers) associated with the site x. On the other hand 7r^(x) denotes 
the number of uuuu dimers, vr^(x) the number of dddd dimers, and itHx) the number of 
oriented (uddu or duud) dimers that live on the bond connecting x and x+u. In our notation 
7T -u{ x ) = 7r u,{ x ~ A) an d similarly for other bond variables. The allowed values for these 
variables are: 

l(x)=0,l n d (x)=0,l n u (x) = 0,l ttJ(x) = 0, 1 tt^) = 0, 1 ^(x) = -1, 0, 1 (2) 

Note that n u (x) = n d (x) = 1 is not allowed since it is absorbed into I(x). Due to the 
Grassmann nature of the fermion fields the following constraints must also be satisfied at 
each site x: 



2I(x) + 



<{x)+<{x)+n u {x)+n d {x) 
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n u {x) + ^ 



x 



nAx) 







(3a) 
(3b) 
(3c) 



where the sum over [i goes over ±1, ... ± (d + 1). Fig.([T]) gives an illustration of an MDPI 
configuration in 1 + 1 dimensions. 

Note that I(x) = 1 on a site x breaks the Ua(X) symmetry but not the SUl{2) x SUr{2) 
symmetry and hence it is called an instanton. In addition to these, each configuration can 
contain neutral pion-loops of alternating uuuu and dddd dimers or charged pion-loops of 
uddu (or duud) dimers. Neutral pions can also form open strings at non-zero m with uu or 
dd monomers at the ends, while charged pions cannot form such open strings in our model. 
Note also that if a bond consists of two dimers (double bond) it can be considered both as 
a neutral pion-loop or a charged pion-loop. Finally, it is useful to recognize that charged 



pion-loops are naturally oriented while neutral pion- loops are not. On the other hand it 
is possible to give an orientation to the neutral pion-loops also. Here we will assume that 
a uuuu dimer on a forward bond from an even site and a dddd dimer on a forward bond 
from an odd site are considered oriented forward. Similarly a uuuu dimer on a forward 
bond from an odd site and a dddd dimer on a forward bond from an even site are considered 
oriented backward. With this convention it is clear that neutral pion-loops are also oriented. 
This will be useful when constructing the algorithm. If we define c = c + m 2 , the partition 
function is given by: 

" " ' (4) 



Z = Y[m nd{x) m nu{x) c I{x) . 

[K] x 



where the sum includes all allowed MDPI configurations. 
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FIG. 1: An example of a 2 x 2 lattice configuration as discussed in the text. 



IV. ALGORITHM 



Since the partition function in Eq. (J3J) is a sum over positive definite weights, a Monte- 
Carlo algorithm can be used to sample the MDPI configurations /C. We have extended the 
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work of [18|, |21[ and constructed the algorithm. It is comprised of four updates : a u <-> d 
/Zip, a /oop swap, and a directed-path mass update and a directed-path fixed-monomer update. 
While each update by itself satisfies detailed balance, one needs a combination of these 
updates to ensure ergodicity. The first three updates are sufficient to ensure ergodicity in 
general. The fixed monomer update in combination with the flip and the swap allows us to 
calculate quantities in a fixed monomer sector which will be useful in some of our results. 
We now discuss each update in detail. 

The u <-»• d flip update This update changes a u quark to a d quark and vice- versa on a pion 
loop or string. Thus, a uuuu dimer becomes a dddd dimer and vice-versa or a uddu dimer 
becomes a duud dimer and vice-versa. Similarly a uu monomer becomes a dd monomer. 
The complete update is as follows: 

1. Let V denote the set of all lattice sites. A lattice site xo e V is selected randomly. 

2. If Xo is part of a pion loop or a string, the u <-> d flip is performed on all dimers and 
monomers associated with the loop or string. The update then ends. 

3. If Xo is part of an instanton or a double dimer, the update ends. 

Because the configuration before and after the update have equal weights, this update au- 
tomatically satisfies detailed balance. 

Loop swap update. This update swaps a neutral pion-loop into a charged pion- loop and 
vice-versa. The complete update is as follows: 

1. A lattice site x G V is selected randomly. 

2. If x is part of a neutral pion-loop, the complete loop is changed into a charged pion- 
loop and vice-versa. The update then ends. 

3. If xo is not part of a loop, the update ends. 

Again, because neutral and charged pion-loops have equal weights, this update automatically 
satisfies detailed balance. 

Directed-path mass update. This update is more complex and can create and destroy 
monomers and instantons while at the same time it can change the shape of the pion loops. 
We construct and use two types of directed-path update that differ on the sites they are 
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allowed to touch. The charged-pion directed path update can only touch sites containing 
either charged pion-loops (including double dimers) and instantons, while the neutral-pion 
directed path update can only touch sites containing neutral pion-loops (including double 
dimers), instantons and monomers. Since the rules of the update are almost the same we 
will discuss them together. We will call the set of sites a particular update is allowed to 
touch as C. 

To understand the directed path update it is useful to think of the MDPI configurations 
as constructed with oriented objects. As we have already discussed above every loop is 
oriented. Open strings of neutral pions are also oriented with an incoming monomer and an 
outgoing monomer at the two ends. An instanton is also an oriented string of zero length 
with a forward and a backward monomer on the same site. A forward monomer is oriented 
into the /i — direction and a backward monomer is oriented from the /i — direction. Thus, 
each MDPI configuration naturally gives every site a forward and a backward direction. The 
directed-path update essentially updates these directions on the sites it touches. 



T 



ru rr 



® 



j — - — ► ® 



FIG. 2: Examples of the active update at a site. The absence of an outgoing arrow implies that 
the update remains on the site. 



The update consists of two parts: an active update and a passive update. During an 
active update one replaces the forward direction ji at the site x by the direction v through 
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FIG. 3: Examples of the passive update at a site. The absence of an ingoing or an outgoing arrow 
implies that the update is on the site. 

which the site was approached. If fi 7^ 0, then after replacing \x with u, the update moves 
to the neighboring site x + £l and replaces the backward direction on the new site with zero. 
Otherwise the update remains on the same site x without disturbing the backward direction 
on it. Some examples of the active update is shown in Fig. (J2J). During a passive update if the 
backward direction // ^ 0, the update ends. Otherwise a new backward direction v is chosen 
with probability w V)X /W. Here w U:X = if x + v C, otherwise w u=0jX = c', w u= ±i = T and 
w V)X = 1 otherwise. W is the sum of weights of all possible choices. One of the v is chosen 
according to heat bath probabilities. When v = is chosen then either a monomer or an 
instanton can be created. If the forward direction on the site is non-zero then the probability 
of creating a monomer is m 2 /c' and creating an instanton is c/c' where d = c + m 2 . If the 
forward direction is zero then only an instanton can be created and so d = c. The update 
ends on a passive site if v = is chosen with a creation of a monomer. On the other hand 



if an instanton is created the update can end only if the forward direction on the site was 
also zero. Otherwise the update remains on that site and continues. If v ^ is chosen then 
the backward direction is changed to v and the update moves to the neighboring site x + v. 
Some examples of the passive update is shown in figure [3J Note that in the charged pion 
directed loop update m 2 is assumed to be zero so monomers can not be generated. 
The complete update is simple and is described below. 

1. A lattice site x G V is selected randomly. If x does not belong to C the whole 
update ends. Further, if both the forward and the backward directions at this site are 
non-zero, the site is active and was reached from the v = direction. Otherwise the 
site is assumed to be passive. 

2. Each current site is updated using the rules of an active site or a passive site succes- 
sively until the update ends. 

Directed path fixed monomer update. This update is identical to the directed path mass 
update, except for a minor change in the neutral pion directed path update. The change is 
such that uu or dd monomers are not allowed to be created or annihilated during the passive 
update. Thus, the update allows the monomers to change positions while keeping the total 
monomer fixed. Note that instanton number can of course change. 



V. OBSERVABLES 



Numerous observables can be measured with our algorithm. The simplest are the three 
helicity moduli or current susceptibilities. In particular, for a conserved current J^{x), the 
helicity modulus (current susceptibility) is defined as: 

X [M=l N X ' ' 

where we are assuming a L t x L d lattice. There are three conserved currents in our model. 
They are the axial, chiral, and vector currents which are given by: 

W = (-l)"K(x) +**{x) + (6a) 

w = (6b) 

W = irUx) (6c) 
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G a Ax,y) = 




ia a (-l)^y) 


G a (x,y) = 






G v (x,y) = 






G a 5 (x,y) = 







The current susceptibilities Y a , Y c and Y v are diagonal observables in our approach and can 
be easily computed for every configuration that contributes to the partition function. 

We can also measure two point correlation functions of many fermion bilinears. In par- 
ticular we can compute 

(7a) 
(7b) 
(7c) 
(7d) 

These correlation functions are non-diagonal observables and cannot be measured easily on 
the configurations that contribute to the path integral. However, they can be measured 
during the charged-pion directed path update as follows: 

1. All the Gi(x,y) are set to zero before the update. 

2. If x is the first active site and y is one of the passive sites visited during the directed 
path update, then Gi(x,y) is changed as follows: 

L L d 

G n (x,y) = G n (x,y) + -^-8 x ^8 VtVO (8a) 

G,(x, y) = G s (x, y) - (-l) x+y ^fs x , Xo 8 y , yo (8b) 

L t L d 

G a {x,y) = G K (x,y)+2yy-6 x>xo 6y >yo (8c) 

G v (x, y) = G 5 (x, y) - (-lf+y^-5 x , xo 5 y , yo (8d) 

If n m (z) is the number of monomers and ni(z) is the number of instantons at any lattice 
site z then the G v and G a correlator get additional "disconnected" contributions: 

L L d 

G a (x, z) = G a (x, z) + [1 + n m (z) + 2n I (z)m 2 /c] ( 9a ) 

G v (x,z) = G v (x,z) - (-l)^[n m (z) + 2n I (z)m 2 /d)}^- (9b) 

(9c) 



3. The function Gi(x, y) thus obtained for each directed path update, when averaged over 
directed path updates, yields the final correlation function defined above. 
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Once the correlation functions are known, the corresponding susceptibilities, Xn, X<?, Xv 
and xs are given by: 



Xi 



x.y 



for i — a, 7T, 77, 5. We have normalized Xi such that 

1 1 d 2 Z 



Xi 



(10) 



(11) 



L d Z (dm) 2 

We have tested our algorithm by comparing the results for all the above observables with 
exact analytic calculations for each of these observables on a 2 x 2 lattice. The exact results 
and those from the algorithm are shown in Appendix I. 



VI. RESULTS IN THE e-REGIME 

In order to establish that we can indeed use our approach to model pions of two flavor 
QCD, in this paper we focus on the e regime of chiral perturbation theory 22 



26 



21, 



23 



24 



28j|. In the phase where chiral symmetry is broken and the anomaly is large, the low 
energy physics of our model must be describable by the Euclidean chiral Lagrangian 

F 2 



C = —trfd^d^ - mT.tr (ll + 



(12) 



where F is the pion decay constant in the chiral limit and E is the chiral condensate and 
U G SU(2) is the pion field. The e regime involves the limit where L, the linear size of the 
four dimensional hypercube, is taken to be large such that FL 1 but mSL 4 is held fixed. 
In this limit a variety of quantities have been computed in the literature. For example the 
behavior of function of L at m = was obtained in [26] for the O(N) model. The 

N = 4 result, relevant here, is 

S 2 L 4 



X* 



3/3i 

1 + 7T77TT7 + 



(FL) 2 (FLY 



+ ... 



The dependence of Y c and Y v were obtained in 27| and it was shown that 



(13) 




+ 0(u 4 



+ 0(u 4 



(14a) 



(14b) 
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for small u = SmL 4 [l + 3/?i/(2(FL) 2 )]. In the above formulas f3\ = 0.14046 is the shape 
coefficient and a,a' b c ,b v are constants that depend on higher order low energy constants of 
the chiral expansion. Also note that Y c = Y v when u = reflects the chiral symmetry of the 
theory. 




20 24 32 40 48 64 

L 



FIG. 4: Vector current susceptibility Yq as a function of L at T = 1.733, c = 0.3 and m = 
[Yq = Yy]. The solid line shows the fit with F = 0.0992 and a' = 2.7. The dotted line shows the 
same curve but with a' = 0. 

We now argue that the calculations in our model are consistent with Eqs. fll3l) and (jT4l) . 
We choose c = 0.3 and T = 1.733 with fixed L t = 2. These parameters are chosen so that 
chiral symmetry is spontaneously broken and the anomalous pion mass {M n ) is about 1 in 
lattice units. As we will see, this choice of T makes F ~ 0.1 in lattice units, which should 
make our results less sensitive to lattice artifacts. In Fig.(j4]) we plot our data for Y c as a 
function of L for m = 0. The solid lines are fits to Eq. (I14|) . The fits are extremely good if 
we allow for a' ^ and use lattice sizes above L ~ 24. We can then extract F = 0.0992(1) 
and a' = 2.7(1) with a x 2 / DOF = 0.8. However, as can be seen from the figure, this means 
that setting a' = 0, i.e.. using only the leading correction in the chiral expansion, will not 
be a good approximation for L < 48. This is clearly due to the smallness of a' / (3\. In 
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other words, although our data is consistent with the Eq. (fT4"]) . unfortunately we are not yet 
sensitive to the universal corrections at order 0(1/ L 2 ). 
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FIG. 5: Chiral condensate susceptibility Xo as a function of lattice size L at T = 1.733, c = 0.3 and 
m = [Xa = Xn]- The solid line is the plot of Eq.(Tj3D with E = 0.1866, F = 0.0992 and a = 3.0. 
The dotted line shows the same curve but with a = 0. 

Let us next focus on the condensate susceptibility Xa- In Fig© we plot x^/L 4 as a 
function of L. The solid line is a fit to the data using the Eq. f|T3|) where we fix F = 0.0992 
and take only the data for L > 20 in the fit. We find S = 0.1862(2), a = 3.0(2) with a 
X 2 /DOF = 1.3. As can be seen again, the universal correction at order 0(1/ L 2 ) is small 
compared to the next order non- universal correction for L < 48. 

In the current work we neglect the log(L) corrections that arise at the order 1/ (FL) 4 
The reason is as follows: Consider for example the chiral condensate susceptibility defined 
Eq. (Hip . Using the results of 3j, |4j, |26j, we can obtain the logarithmic corrections to Eq. (1131) . 
In particular we find that 



S 2 L 4 



Xn 



3ft 
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where now a = (3/3? + 15/3 2 )/2 + 3[log(A M /F) + 41og(A s /F)]/167r 2 and (3 2 = -0.020305 is 
another shape coefficient. The mass scales Am, As encode the non-universal information of 
our model and are defined in 26j. Assuming L\ = 20 is the smallest lattice size and L 2 = 64 
is the largest lattice size we use in the fits, we see that the change in the logarithmic correction 
term 151og(L 2 /-Li)/(167r 2 ) ~ 0.1 is within errors of the constant a = 3.0(2) obtained above 
by fitting the Xn data to Eq. (fT3~]) . Thus, we believe our errors are still large and we are 
not yet sensitive to the logarithmic corrections. Interestingly, since 15 log(FL)/(167r 2 ) is 
much smaller than a in the region of our fits, we can estimate that a ~ 3[log(Ajvf/F) + 
41og(A E /F)]/(167r 2 ), which means that [log(A M /F) + 41og(A s /F)] (~ 150) is unnaturally 
large, and the factor is essential to keep the coefficient of 1/(FL) A of order 1. In other 
words factors like —r^ cannot always be assumed to be small since they can be multiplied 
by large numbers. 

Having confirmed that our results are consistent with Eq. (jT4")) when u = 0, we can also 
check that our model gives results consistent with Eq.( fT4l) at order u 2 . One way to do 
this is to tune the quark mass and the volume such that u is fixed and small. This is 
cumbersome. For example even at u = 1 we see that for L = 48 our quark mass should be 
as small as 10~ 6 . Since m 2 is involved in a probability, one begins to worry about double 
precision arithmetic. Thus, here we devise another method. To understand our approach 
let us expand the partition function in powers of the quark mass 

Z = Z + m 2 Z 2 + m 4 Z 4 + .... (16) 

where Z n is obtained from configurations with n monomers. In this expansion we can 
neglect the m 2 contribution to instanton weights as they will not contribute in the e regime. 
Similarly, let Yy^ and be the value of the current susceptibilities when computed in 
the n monomer sector. Expanding observables in the various monomer sectors it is possible 
to show 

Y5 = yf +m'{yf>-yf }§ + ... (it) 

where i = v,c. From Eq. (|T3l we see that 



f- s x( 1 + w + (^ + -) (18) 
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Since at u = no monomers contribute we also have 



Substituting our knowledge of Z 2 /Z and F c ° we get 



«=IV<i + 74 i + 4 + £l*f-n m > + -l ™ 



(FL) 

Comparing with Eq. pijl we conclude that 

y — t( 1 + W + '') (21b > 

(2) (2) 

Notice that Yc ^ Yv due to the effects of explicit chiral symmetry breaking that is 
introduced due to the presence of the monomers. 




(2) (2) 

FIG. 6: Plot of Y C K ' and Yj J , evaluated in the two monomer sector as a function of L at T = 1.733, 



c = 0.3 and m = 0. The solid lines are fits to Eq. ([21 j) discussed in the text. 
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We have used the fixed monomer update to compute Y} and Y} . In Fig. ([6]) we plot 
our results as a function of L. We fix F = 0.0992 and try to fit our data to Eq. (l2T|) . We 
obtain b' c = 4.1(1) with a x 2 /DOF = 1.1, and b' v = 4.2(1) with a x 2 /DOF = 2.1. Although 
our results again appear consistent with the predictions at large L, the large values for the 
constants b' c and b' v show that we need data with small errors for L > 48 to be sure we can 
be sensitive to the universal predictions at 0(1/ L 2 ). 

While our model is analogous to real QCD, we can not expect our results to be equivalent 
to real QCD. However, for pedagogical reasons, we can make comparisons to real QCD 
studies. In particular, at T = 1.733, we found Fa\ at = 0.0992 and M v a\ at = 1.1 where we 
have introduced ai at to be the lattice spacing. If we take ai at = IGeV -1 , we find F ~ 99MeV 
and M„ ~ l.lGeV which are reasonably close to physical values. 




FIG. 7: F and £ plotted as a function of T. The solid lines are fits to Eq. (|22p with Ap = 0.943, 
A s = 1.769 and T c = 1.73779. 

VII. CONTINUUM LIMIT 

At c = 0.3 and T = 1.733 we found that F ~ 0.1. Such a small value of F clearly is 
not an accident. At a fixed c ^ the model should go through a phase transition at some 
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T 


F 


a X 2 


S 


b X 2 


M v 


x 2 


1.68 


0.2711(1) 


1.28(3) 0.3 


0.5045(2) 


1.22(4) 1.1 


1.04(5) 


1.5 


1.70 


0.2327(1) 


1.91(9) 1.2 


0.4346(2) 


1.72(11) 0.2 


1.06(5) 


1.2 


1.72 


0.1744(1) 


2.58(13) 0.1 


0.3268(2) 


2.44(17) 0.5 


1.02(3) 


1.1 


1.73 


0.1234(2) 


2.25(10) 1.7 


0.2318(2) 


2.26(8) 1.5 


1.05(2) 


1.6 


1.7325 


0.1038(1) 


2.80(18) 0.1 


0.1947(2) 


2.77(23) 0.7 


1.08(3) 


1.3 


1.733 


0.0992(1) 


2.71(9) 0.8 


0.1866(4) 


3.00(20) 1.8 


1.05(3) 


1.5 


1.734 


0.0889(3) 


2.78(15) 2.3 


0.1672(4) 


2.81(17) 2.3 


1.08(3) 


1.4 



TABLE I: Fits of Y c and x<j as a function of L based on Eqs. (|14p and (|13p at c = 0.3 and m = 
at various values of T < T c . The last two columns give the value M v and the x 2 /DOF obtained 
from a fit to G v {x,y). 



critical temperature T c . If this transition is second order then at the critical temperature 
one would expect F to vanish. Thus, T = 1.733 must be close to such a critical point. Since 
we are in four dimensions, we expect quantities to obey mean field scaling near a second 



order critical point. In particular this means 29]: 



F ~ A F (T c - T) -2 1 log(T c - T) 1 3 (22a) 
£ ~ A S (T C - T) l 2\log(T c - T)|3 (22b) 

We have verified that our results for c = 0.3 are consistent with this expectation. Tab.(JI]) 
gives the chiral perturbation theory fit parameters as a function of T. The values of F and £ 
are plotted as a function of T in Fig. (JTj). If we use the data for T > 1.73 a combined fit of F 
and £ as a function of T to Eq.O gives A F = 0.943(4), A s = 1.769(4) and T c = 1.73779(4) 
with a x 2 /DOF = 0.7. 

Although we cannot rule out a weak first order transition with a large correlation length, 
if we assume that there is a second order critical point at T = T c , we can in principle make F 
arbitrarily small by tuning our T closer and closer to T c . This raises an interesting question: 
By approaching T c can we define a "continuum limit" of our lattice pion theory by holding 
F to be the fixed physical scale? However, due to trivi ality of four dimensional scalar field 



theories this limiting theory is expected to be free 30j, |3l|, [32| . On the other hand, chiral 



perturbation theory suggests that pions will always interact as their momentum increases 
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which is in contradiction to the statement of triviality. Hence, we conjecture that the mass of 
the a resonance will go to zero in units of F. Since triviality is a logarithmic effect we expect 
M a jF ~ 1/| log(T c — T)]\ p where p is some power. This means that the chiral expansion 
would break down at momenta close to M a rather than F and the usual power counting 
of chiral perturbation theory would become questionable close to T c . Thus, our approach 
to model pion physics of QCD should eventually fail! On the other hand since triviality is 
a logarithmic effect, it is likely that we still will have a window where our model produces 
pions very much like in two-flavor QCD. The present study is an attempt to demonstrate 
this in the e-regime. 



VIII. CONCLUSIONS AND FUTURE WORK 



In this work we have developed a new approach to model the physics of pions in Nf = 2 
QCD. Our approach uses Nf = 2 lattice QED in the strong coupling limit. We have shown 
that using the mapping to dimer models we can study our model very efficiently in the 
chiral limit and close to it. We have established consistency with the e regime predictions of 
chiral perturbation theory. We have also demonstrated that we can make F <C 1 by tuning 
a fictitious temperature so one approaches a second order phase transition. This tuning 
is expected to help remove lattice artifacts and approach a continuum-like theory. But, 
as explained in section IVHl our method will eventually fail due to logarithmic triviality of 
scalar field theories. Since triviality is a logarithmic effect, there should be a large window 
where our method may provide a model for pions. Our results confirm this belief. 

There are many directions for the future. We are currently performing calculations in 
the p-regime and hope to understand the re gion in our mode 
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35, 



36 



37 



where chiral perturbation 
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,[3J 



Since we have 



theory, up to a particular order, will be valid [33j, 
already computed the leading low energy constants in the e-regime we can constrain them 
when we analyze the p-regime. Another interesting arena to explore is pion scattering and 
resonance physics by measuring the appropriate two and four point correlation functions and 



extracting scattering lengths and phase shifts via Liischer's method 
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421 ] . Our approach 



should also allow one to study the effects of the quark mass on these quantities. 
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Appendix I 



Here we give the analytic expressions for the partition function, the chiral, vector, and 
axial helicity moduli and the two susceptibilities for a 2 x 2 lattice, as defined in SectionfTVT). 

(Note c = c + m 2 ): 

Z(T, c, m) = 36T 4 + 64T 2 + 36 + c 4 + 12(1 + T 2 )c 2 + 8(1 + T)c 2 m 2 (23a) 

+ 32Tcm 2 + 16(1 + T 2 )m 4 + 48(1 + T 3 )m 2 + 32(T + T 2 )m 2 

2Z x Y$(T, c, m) = 96 + 64T 2 + 16c 2 + 112m 2 + 32m 4 + 8c 2 m 2 + 32Tcm 2 (23b) 

+ 32(2T + T 2 )m 2 

2Z x y£(T, c, m) = 96 + 64T 2 + 16c 2 + 64m 2 (23c) 

Y*(T, c, m) = Y°(T, c, m) + (23d) 

2Z xxn = 24T 3 + 16T 2 + 16T + 24 + 4(1 + T)c 2 + 4(3 + 4T + 3T 2 )c + 2(S3e) 

+ 8(1 + T)cm 2 + 16(1 + T + T 2 )m 2 

2Z x Xr, = 24T 3 + 16T 2 + 16T + 24 + 4(1 + T + m 2 )c 2 - 4(3 - 4T + 3T 2 |23f) 

- 2c 3 - 8(1 + T)cm 2 + 8(5 + 6T + 5T 2 )m 2 + 16(1 + T)m 4 



In Tabs.f UIj) and fllllj) we compare the analytic results with the results obtained using our 
algorithm for different parameters. The agreement gives us confidence that our algorithm 
must be correct. 
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T c m Algo. Exact Algo. Exact Algo. Exact 











Y 

± w 






w 






Y v 

w 




1.0 


0.5 








0.8023(9) 





80246.. 


. 0.5763(6) 





57721.. 


. 0.5771(8) 


0.57721... 


1.5 


0.5 








0.5212(7) 





52141.. 


. 0.3275(6) 





32790.. 


. 0.3274(7) 


0.32790... 


1.0 


1.0 








0.7449(8) 





74534.. 


. 0.5470(7) 





54658.. 


. 0.5468(8) 


0.54658.. 


1.5 


1.0 








0.4967(7) 





49720.. 


. 0.3178(5) 





31821.. 


. 0.3175(6) 


0.31821... 


1.0 


1.5 








0.6667(8) 





66645.. 


. 0.5016(5) 





50240.. 


. 0.5014(6) 


0.50240... 


1.5 


1.5 








0.4607(6) 





46147.. 


. 0.3029(5) 





30325.. 


. 0.3033(6) 


0.30325... 


1.0 


2.0 








0.5812(7) 





58064.. 


. 0.4519(6) 





45161.. 


. 0.4521(6) 


0.45161... 


1.5 


2.0 








0.4199(5) 





41927.. 


. 0.2853(5) 





28451.. 


. 0.2853(5) 


0.28451... 


1.0 


0.0 





5 


0.7734(5) 





77336.. 


. 0.5972(5) 





59730.. 


. 0.4867(6) 


0.48692... 


1.0 


0.0 





r 




U.5Uoz(6 J 


u 


5U lU2.. 


. U.34 ly(o) 


u 


34034.. 


. U.2oz/(4J 


U.Z0319... 


1.25 


0.5 


1 


5 


0.4192(4) 





41948.. 


. 0.3914(6) 





39164.. 


. 0.1285(3) 


0.12831... 


1.5 


1.25 


5 


0.4516(5) 





45126.. 


. 0.3250(4) 





32496.. 


. 0.2564(4) 


0.25611... 


2.0 


1.25 1 


5 


0.2664(3) 





26636.. 


. 0.2395(3) 





23928.. 


. 0.0733(1) 


0.07318... 


1.0 


0.5 


2 





0.3834(5) 





38411.. 


. 0.3757(4) 





37622.. 


. 0.0915(2) 


0.09122 



TABLE II: Helicity moduli for a 2 x 2 lattice as discussed in the text. 
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XlT Xr] 

T c m Algo. Exact Algo. Exact 



1.0 


0.5 











3601(3) 


0.359877.. 


. 


2172(2) 


0.217334... 


1.5 


0.5 











2676(2) 


0.267764.. 


. 


1824(1) 


0.182429... 


1.0 


1.0 











4040(3) 


0.403727.. 


. 


1429(2) 


0.142857... 


1.5 


1.0 











2981(2) 


0.298322.. 


. 


1367(1) 


0.136731... 


1.0 


1.5 











4220(3) 


0.422301.. 


. 


0801(1) 


0.080103... 


1.5 


1.5 











3172(3) 


0.317264.. 


. 


0948(1) 


0.094767... 


1.0 


2.0 











4194(2) 


0.419355.. 


. 


0323(1) 


0.032258... 


1.5 


2.0 











3251(3) 


0.324754.. 


. 


0588(1) 


0.058961... 


1.0 


0.0 





5 





2849(2) 


0.28481... 





2849(2) 


0.28481... 


1.5 


0.0 





5 





2220(1) 


0.22220... 





2220(1) 


0.22220... 


1.25 


0.5 


1 


5 





1708(2) 


0.17081... 





1412(2) 


0.14121... 


1.5 


1.25 


5 





2767(2) 


0.27650... 





1175(1) 


0.11753... 


2.0 


1.25 1 


5 





1384(1) 


0.13831... 





0933(1) 


0.09323... 


1.0 


0.5 


2 








1339(1) 


0.13365... 





1160(1) 


0.11584... 



TABLE III: Susceptibilities for a 2 x 2 lattice as discussed in the text. 
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